Pathlength dependence of energy loss within in-medium showers 
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Studying the pathlength dependence of high Pt hadron suppression in heavy-ion collisions by 
measuring the dependence of hard hadron production on the angle <j) with the reaction plane in non- 
central collisions has so far been one of the most successful tools in constraining the microscopical 
picture of leading parton energy loss. With the imminent start of the LHC heavy-ion program 
and the possibility of full jet reconstruction, the focus has shifted to models which are capable of 
simulating full in-medium parton showers rather than tracing the leading parton only. Yet, on the 
level of single inclusive hadron observables, such shower models need to reproduce the findings of 
leading parton energy loss models, in particular they need to reproduce the correct reaction plane 
angle dependence of hadron suppression. The aim of this work is to illustrate at the example of the 
in-medium shower code YaJEM (Yet another Jet Energy-loss Model) how pathlength dependence 
arises in a MC shower formulation, how it relates to pathlength dependence of leading parton 
energy loss and how model results compare with data on the reaction plane dependence of the 
nuclear suppression factor Raa^)- 

PACS numbers: 25.75.-q,25.75.Gz 
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I. INTRODUCTION 

The idea to use the modification of high Pt hadron 
production due to the bulk medium produced in ultra- 
relativistic heavy-ion collisions to study properties of the 
medium is known as 'jet tomography' [lHa]. However, 
in order to gain any tomographical information about 
the macroscopical medium density distribution and its 
evolution from measured data in practice, it is necessary 
to know first the microscopical mechanism by which the 
medium modifies high Pt hadron spectra. A less ambi- 
tious approach is to determine under what assumptions 
a combination of macroscopical evolution in terms of a 
bulk evolution model like fluid dynamics and microscopi- 
cal parton-medium interaction model can describe a large 
body of data in a systematic way So far, most of 

these studies have been carried out for models which de- 
scribe the medium modification as energy loss from the 
leading shower parton prior to hadronization. 

An emerging result is that the pathlength dependence 
of energy loss, accessible experimentally e.g. via the de- 
pendence of the nuclear suppression factor Raa{4>) on 
angle of outgoing hadrons with the reaction plane <f> pro- 
vides a crucial constraint for the microscopical physics. 
In particular, the pathlength dependence cannot be lin- 
ear in a constant medium (as characteristic for elastic or 
generally incoherent processes) [IJ [H[ but must be at 
least quadratic II 311 characteristic of coherence, possibly 
even cubic fl3l . Il4l | as found in some versions of strongly 
coupled QCD-likc models. 

On the other hand, going beyond the capability to 
describe single inclusive hadron observables to full jets, 
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Monte-Carlo (MC) models capable of describing the de- 
velopment of a parton shower originating from a highly 
virtual initial parton inside a medium, followed by 
hadronization into jets are being developed. Both elastic 
collisions [To| as well as medium-induced radiation 

have been considered as mechanisms for the modifi- 
cation. 



Monte-Carlo (MC) models are not natural frameworks 
to incorporate the physics of quantum coherence (which 
appears necessary to reproduce the pathlength depen- 
dence of energy loss seen in the data), but it has been 
demonstrated in [2fj| that an algorithm can be devised 
which consistently reproduces known eikonal infinite en- 
ergy analytical limits. Yet in the same work, it has also 
been shown that including realistic effects such as finite 
energy corrections, the pathlength dependence is effec- 
tively weaker than quadratic (see also [HI). This is in 
manifest contradiction to systematic investigations which 
show that the data for (</>) are marginally in agree- 
ment with a quadratic dependence and would be easier 
reproduced with a cubic dependence [lH . 



The aim of this paper is to illustrate the connection be- 
tween the pathlength dependence of energy loss in lead- 
ing parton models and the effect of pathlength on an 
in-medium shower more clearly and to suggest a possi- 
ble resolution for the apparent puzzle outlined above by 
suggesting that the pathlength dependence observed in 
the data might be due to a yet different effect. For all 
computations in the following, the radiative energy loss 
scenario (RAD) of the MC in-medium shower code Ya- 
JEM is used. 
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II. THE MODEL 



A. General description 



medium is increased by the ratio of their Casimir color 
factors, 3/ 1 = 2.25. The medium is characterized via the 
parameter q(() where a relation like 



YaJEM is described in detail in [ly, [l7| , here we just 
reproduce the essential details necessary to understand 
how pathlength dependence arises within the model. 

In the absence of a medium, YaJEM utilizes the 
PHYSHOW algorithm [22j in momentum space to model 
the evolution from a highly virtual initial parton to a 
shower of low virtuality partons and the Lund model [23 1 
to hadronize this shower. In this framework, the QCD 
evolution from a highly virtual state to a parton shower 
is modelled as a series of branchings a — > be where the 
virtuality scale decreases at each step. 

In the presence of a medium, it is no longer sufficient 
to model in momentum space only since the medium evo- 
lution is caused by pressure gradients in spacetime and 
the density of matter felt by hard partons is a function 
of position space rather than momentum space. In order 
to make a link between the shower evolution in momen- 
tum space and the medium density evolution in position 
space, YaJEM utilizes the uncertainty relation to com- 
pute the average formation time of a virtual daughter 
parton b generated in the branching of a parent a as 



(i) 



where Ei denote parton energies and Qi parton virtu- 
alities. YaJEM provides the additional option not to use 
the average formation time but rather for each branching 
an actual formation time given by an exponential distri- 
bution 



P(n) = exp 



n 

(n) 



(2) 



which averages to (77,). This allows to specify the tem- 
poral structure of the shower, and combined with the 
assumption that all shower partons sample the medium 
along the eikonal trajectory of the shower-initiating par- 
ton, the complete space-time picture of the shower is 
fixed. 

In the radiative energy loss scenario (RAD) of YaJEM, 
the parton-medium interaction is modelled by a medium- 
induced virtuality gain during the parton formation time 



q(Q =K-2- [e(C)] 3/4 (coshp(C) - sinhp(C)cosV) (4) 

is assumed as the link to fluid dynamical parameters 
such as the energy density e or the local transverse flow 
rapidity p (at angle with the parton trajectory). 

One will observe that Eq. (lU),© arc not indepen- 
dent — parton virtuality determines the parton forma- 
tion time, but the formation time in turn determines the 
virtuality gain from the medium. If AQ 2 a <C Q 2 , holds, 
i.e. the virtuality picked up from the medium is a cor- 
rection to the initial parton virtuality, we may add AQ 2 
to the virtuality of parton a before determining the kine- 
matics of the next branching. If the condition is not ful- 
filled, the lifetime is determined by Q 2 + AQ 2 and may 
be significantly shortened by virtuality picked up from 
the medium. In this case we iterate Eqs. CD),© t° deter- 
mine a self-consistent pair of ((r a ),AQ^). This ensures 
that on the level of averages, the lifetime is treated con- 
sistently with the virtuality picked up from the medium. 
The actual lifetime is still determined by Eq. ([2]). This 
consistency condition is fundamentally the same as the 
algorithm presented in [20]. It corresponds to the limit 
of a coherent summation of a large number of medium 
interactions during the formation time (in YaJEM, virtu- 
ality gain is a smooth function of pathlength, rather than 
a sum over discrete interactions transferring virtuality) 
and hence represents LPM interference in this limit. 

However, for technical reasons, despite the consistency 
condition for medium-induced virtuality gain, the in- 
medium shower simulated by YaJEM remains an ex- 
pansion for thin media. The reason is that, as in the 
PYSHOW algorithm, parton splitting no longer occurs 
once a parton has reached a lower virtuality scale Qo 
(where Q = 0(1) GcV — to be more precise, for mas- 



sive partons the limit is given by m q e f f 
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with 
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iQo and m q 



eff 



q + jQo anC ^ m q tne CUr " 

rent quark mass). Partons falling at any point of the 
evolution below this scale are set on-shell and do not 
longer branch. YaJEM inherits this condition also in a 
medium. This limits the distance scale L max for which 
an in-medium shower is evolved reliably paramctrically 
to L max ~ E/Qq where E is the energy of the shower- 
initiating parton. 



AQ* = / dCq{C) (3) 

Jt° 

where the time r a is given by Eq. the time t° 
is known in the simulation as the endpoint of the pre- 
vious branching process and the integration d£ is along 
the eikonal trajectory of the shower-initiating parton. If 
the parton is a gluon, the virtuality transfer from the 



B. Pathlength dependence 

It should now be clear from where any dependence of 
parton distributions created by YaJEM on the in-medium 
pathlength arise. 

• The leading effect is given by Eq. © which states 
that the amount of medium-induced virtuality is 
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a monotonically increasing function of the path- 
length. If q is a constant, this states that the to- 
tal amount of medium-induced virtuality (or the 
number of scattering centers probed by a parton) 
is linearly proportional to the pathlength. Since 
Eq. ([3]) holds for every parton in the shower, this 
would lead to a crescendo-like scenario [24| for the 
virtuality transfer to the shower as a whole in the 
absence of finite energy or finite virtuality correc- 
tions (in practice however, such corrections turn 
out to dominate the dynamics (25j). 

• The consistency condition between formation time 
and medium-induced virtuality forms a second im- 
portant ingredient. It states that low Q 2 , i.e. 
collinear medium-induced radiation may not occur 
after arbitrarily short formation times. When ap- 
plied to a single almost-on-shell parton propagat- 
ing through a constant medium, it gives rise to a 
quadratic pathlength dependence of mean radiated 
energy as in [2(3]. However, note that since YaJEM 
requires Q > Qo, this situation is strictly never re- 
alized. To some degree, its effect is also obscured 
by the randomization of formation times Eq. ([5]) as 
will be demonstrated later. 

• The minimum virtuality scale Qo down to which 
partons are evolved induces a spurious length de- 
pendence whenever a pathlength L > L max is 
probed — for L 3> L max (i.e. far outside the valid- 
ity region of the model), all partons in the shower 
are on-shell and no branching occurs for any den- 
sity of the medium. Thus, when applied outside the 
region of validity, YaJEM leads to less modification 
of showers than expected in reality 

• An interesting proposal, suggested in (2(| |27| , is to 
define the minimum virtuality scale by 



Qo = y/E/L (5) 

where L is the actual pathlength a shower is evolved 
in the medium and E the energy of the shower- 
initiating parton, rather than leave it a constant. 
If this is done, it replaces the spurious Qo depen- 
dence by a genuine scale dependence of the shower 
evolution which states that the medium cannot in- 
duce radiation below a scale Qo- We will test this 
prescription in the following. 

In practical terms, it was found in [l6l [l7j that many 
observablcs (e.g. the fragmentation function) follow to 
good approximation a scaling law with the integrated to- 
tal virtuality gain from the medium, AQ 2 = J Q d^q(f) 
inside the region of validity. This implies that the path- 
length dependence of such observables is approximately 
linear in a constant medium. 



C. Implications of a dynamical minimum virtuality 

Let us take some time to consider the implications of 
Eq. ([5]) in some more detail. The obvious problem is that 
Qo regulates both the endpoint for the vacuum and for 
the in-medium parton shower evolution (which can't be 
disentangled on a branching- by-branching basis). While 
length dependence in the medium contribution is reason- 
able, there is no justification for any length effect on the 
vacuum shower evolution. In contrast, in vacuum, Qo 
regulates the transition between the perturbative par- 
tonic and non-pcrturbative hadronic domain of fragmen- 
tation. This in turn is an ill-defined condition in a de- 
confined medium, as hadrons cannot be formed in such 
a medium. 

Eq. ([5]) can therefore not be useful in general but must 
be seen as an approximation for a situation in which the 
shower evolution is dominated by medium effects and 
in which the kinematics is such that the possibility of 
hadrons being formed in the medium is not a crucial is- 
sue. The latter condition is readily fulfilled when consid- 
ering single inclusive high Pt pion observables (as done 
in this work), as the pion formation length (estimated 
as t ~ E n /m^) is far larger than the medium size. Note 
that this is not the case for the low z part of the fragmen- 
tation function or for the production of heavier hadrons. 




z 



FIG. 1: (Color online) Response of the computed fragmen- 
tation function of a 20 GeV quark into charged hadrons to 
choices of the minimum virtuality scale Q% from 0.5 GeV 2 to 
2 GeV 2 , both for the vacuum case and for a medium-modified 
shower for total medium induced firtuality AQ 2 = 6 GeV 2 as 
typical for the majority of partons in RHIC high Pt observ- 
ables within YaJEM. 

It turns out that the first condition is also fulfilled for 
typical in-medium pathlengths as relevant for computing 
observables. FiglT] shows vacuum and medium-modified 
fragmentation functions as bands indicating their re- 
sponse to a change in Qq from 0.5 GeV 2 to 2 GeV 2 for 
typical parton paths at RHIC kinematics. It is readily ap- 
parent from the figure that the majority of the variation 
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with Qq seen in the in-medium case cannot be accounted 
for by the comparatively small effect of the scale change 
in vacuum, indicating that the pathlcngth condition in 
this situation indeed dominates over the scale dependence 
of the transition to the non-perturbative regime. Part of 
the reason is that Qo does not only define the endpoint of 
the perturbative shower description, but also the starting 
point of the non-perturbative Lund hadronization model. 
Thus, to some degree, in the vacuum decreased showering 
is compensated by an increased action of the Lund model, 
whereas no such compensation occurs in the medium case 
(since the hadronization is by assumption always treated 
as in vacuum). 

Nevertheless, Eq. ([5]) has to be regarded as a first ap- 
proximation. Ultimately, a criterion which also takes into 
account the virtuality AQ 2 induced by the medium along 
the parton path would be superior. This can be under- 
stood in the limit of small pathlcngth L. For applications 
to observables, AQ 2 and L are not independent and small 
L also implies small AQ 2 . In this case, the appropriate 
limit for Qq is the vacuum default value Qo = 1 GeV 
rather than the value given by Eq. ([5]), as the assump- 
tion that the medium effect dominates the dynamics of 
the shower is not fulfilled when AQ 2 is small. However, 
one could think of an artificial problem in which AQ 2 is 
large while L is small, in which case the limit for small 
pathlcngths is not the vacuum. 

For the time being we will nevertheless explore the 
consequences of the approximation Eq. (JS|) as it stands 
rather than introduce a more general prescription which 
brings additional model dependence. This means that 
the results are not reliable for small pathlengths where 
the shower is not dominated by the medium, but the un- 
derlying assumption is justified by the vast majority of 
partons encountered in the computation of experimen- 
tally relevant observables. 



D. Medium-modified showers and energy loss 

It is not straightforward for medium-modified shower 
models to make contact with leading parton energy loss 
models. Typically in an energy loss model like [l|-|3| there 
is a clear distinction between the leading parton and radi- 
ated partons made by the a priori assumtion that the en- 
ergy of the leading parton is much larger than the energy 
of radiation quanta. In this limit, the medium-induced 
energy loss is given by the sum of the energy of all ra- 
diated quanta or equally by the energy difference of the 
leading parton before and after the passage through the 
medium. A posteriori, one may however find (rare) pro- 
cesses in which the energy of a radiated parton is much 
larger than the energy of the parent — these are formally 
counted as large energy loss. 

In a shower model, the leading (most energetic) parton 
can be identified at all times, but the identity of that 
parton may change easily during the evolution. Unlike 
in the energy loss case, the radiation of a hard gluon 



carrying away 80% of the energy of a parent quark would 
not be counted as large energy loss but as change in the 
identity of the leading parton. No particular branching 
process can be attributed uniquely to vacuum or medium 
— differences between a medium-modified and a vacuum 
shower can only be defined at the level of averages over 
many events. This implies that the baseline to compare 
with is not a single parton with a definite energy, but 
rather a vacuum shower with a statistical distribution of 
hadrons at given momenta. 

As a way to overcome these difficulties and to make 
contact between in-medium shower and energy loss mod- 
elling, in [T3| c-quarks as shower initiators have been con- 
sidered. Due to the deadcone effect, these have compar- 
atively low vacuum shower activity and a hard fragmen- 
tation pattern, thus minimizing the distortion of the vac- 
uum baseline from a single parton with definite energy. 
Moreover, since c quark generation in branchings g — > cc 
is mass suppressed, tagging the c quark almost always 
tags the shower initiator. 

Energy loss of the leading c quark can then be deduced 
from the ansatz 

j at med /> rl AT vac 

— c (E) = J d(AE)—jE')P(AE)6(E'-E-AE) 

(6) 

solving for the probability distribution of energy loss 
P(AE). Note that this ansatz contains the rather dras- 
tic assumption that there is no parametric dependence 
on the initial energy E. If we require P(AE) to be a 
probability distribution, the assumption may imply that 
for some partons in the distribution dN/dE^ ac the en- 
ergy loss AE is larger than their energy E in which case 
they have to be considered lost to the medium. 

In practice, Eq. ((6j) is discretized and cast it into the 
form of a matrix equation 



N^E?) = K l3 (E\ AE^iAW) (7) 

3=1 

where dN/dE c is provided at to discrete values of E 
labelled Ni and P(AE) is probed at n discrete values of 
AE labelled Pj. The kernel Kj is then the calculated 
dN/dE c for all pairs (E i ,AE : >) where the energy loss 
acts as a shift of the distribution, i.e. dN/dE™ ed (E) = 
dN/dEl ac {E + AE). Eq. (JTJ) can in principle be solved 
for the vector Pj by inversion of Kij for to = n. How- 
ever, in general this does not guarantee that the result 
is a probability distribution. Especially in the face of 
statistical errors and finite numerical accuracy the direct 
matrix inversion may permit negative Pj which have no 
probabilistic interpretation. Thus, a more promising so- 
lution which avoids the above problems is to let to > n 
and find the vector P which minimizes ||iV — J\P|| 2 sub- 
ject to the constraints < Pi < 1 and Y^i=\Pi = 1- 
This guarantees that the outcome can be interpreted as 
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a probability distribution and since the system of equa- 
tions is overdetermincd for m > n errors on individual 
points Ri do not have a critical influence on the outcome 
any more. 

Note that this ansatz leads to additional complications 
in combination with the prescription to set the minimum 
virtuality scale Qq in the shower Eq. ([5]). This is due 
to two reasons: First, in essence for short paths, energy 
loss cannot be extracted reliably because the scale choice 
for the vacuum part becomes only irrelevant if the vac- 
uum shower is a small correction to the medium-induced 
shower as discussed above. Unfortunately, for truly long 
paths the second complication becomes apparent: unlike 
in the light quark case, the charm mass m c rather than 
Eq. (|5|) determines the point at which the branching is 
stopped. In order to avoid the second limit, we use very 
energetic c quarks to extract energy loss, thus allowing 
for a larger length scale to be probed. 

For these reasons, the actual effect of Eq. §5§ on the 
pathlength dependence of light parton energy loss is 
only in a very limited way represented by evaluating 
Eq. (J7J, and all energy loss results shown for this pre- 
scription need to be treated with the appropriate amount 
of caution. However, this is largely a complication of 
the way energy loss is determined from an in-medium 
shower result — there is no compelling reason to suspect 
that Eq. ([5]) would be as problematic for the in-medium 
shower evolution result as such or in the computation of 
observables. 

As a side remark, the fact that the mass of heavy 
quarks largely set the scale at which the shower evolu- 
tion is terminated implies that there is no substantial 
heavy-quark radiative energy loss in the kinematical re- 
gion of RHIC experiments in YaJEM. Without an ad- 
ditional mechanism like elastic energy loss, YaJEM thus 
fails (like other radiative energy loss models) to account 
for the suppression of heavy quarks seen in single electron 
spectra. 



III. RESULTS 

A. Pathlength dependence of energy loss 

In Fig. [5] the mean energy loss as a function of path- 
length extracted from YaJEM via Eq. for a 100 GeV 
charm quark as shower initiator is shown for different 
options. 

In the first case, the formation time of virtual states 
Eq. (TT|) is not randomized by Eq. @. The result is ini- 
tially very well described by a quadratic length depen- 
dence until (around 3 fm) finite energy corrections be- 
come important and the limit of applicability is reached 
and turn the curve over. 

If randomization is taken into account the quadratic 
dependence becomes much blurred, and the net result 
can be described well by a linear growth with pathlength 
until finite energy corrections and the limit of applicabil- 



100 GeV c quark, q = 5 GeV 2 /fm 




L [fm] 

FIG. 2: (Color online) The integrated mean energy loss AE 
as a function of pathlength extracted from YaJEM for a 100 
GeV shower-initiating c quark propagating through a con- 
stant medium with q = 5 GeV 2 /fm for different options 
(see text), compared with a quadratical dependence on path- 
length. 



ity turn the result over to a constant again. This agrees 
nicely with the scaling law identified in [l~7j and shows 
that an algorithm which in principle reproduces the an- 
alytical L 2 dependence in the right limit may not show 
this behaviour in practically relevant situations. 

Finally, the last case shows the effect of a variable min- 
imum virtuality scale Qq according to Eq. ([5]). Note that, 
following the previous discussion, certainly beyond 10 fm 
the result is not characteristic for what would happen 
in the light quark case as the charm quark mass rather 
than Qq determines where the branching process is ter- 
minated. For a light quark, the curve would rise more 
strongly. 



B. Pathlength dependence of the nuclear 
suppression factor Raa 

In order to illustrate how pathlength effects manifest 
in observable quantities, we use a realistic 3+ld hydro- 
dynamical evolution model [28| for the bulk medium and 
compute the dependence of the nuclear suppression fac- 
tor Raa on the angle with the reaction plane in this 
background. The strai ghtf orward computation is de- 
scribed in detail in pH l29j|. It involves averaging the 
medium-modified fragmentation function (MMFF) com- 
puted from YaJEM over all possible paths through the 
medium, which are in turn determined by the distribu- 
tion of binary collisions in the transverse plane and the 
orientation of the parton momentum vector with respect 
to the reaction plane. This result is then folded with 
a leading order pcrturbative Quantum Chromodynam- 
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ics (LO pQCD) calculation for the parton spectrum to 
get the yield of high Pt hadrons. Since YaJEM obtains 
the MMFF at given partonic scale whereas LO pQCD 
computations require the fragmentation function to be 
known at given hadronic scale, a scale matching proce- 
dure as described in [2j| has been used in the following. 

In all computations, the factor K in Eq. (|4|) is adjusted 
such that the measured Raa f°r 0-10% central Au-Au 
collisions at P T = 6 GeV is reproduced by the calcula- 
tion. Both the Pt dependence of Raa, the centrality 
dependence and the dependence on the angle <f> with the 
reaction plane follow then without additional free param- 
eters. Note that the value of K depends on the choice of 
Qo (as pointed out already in [I?} ) ■ 
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FIG. 3: (Color online) The nuclear suppression factor Raa 
for 0-10% central 200 AGeV Au-Au collisions computed for 
two different options (see text) in YaJEM, compared with 
PHENIX data 

In Fig. [3] the result for the angular average of Raa is 
shown in comparison with PHENIX data [30] for both a 
fixed value of the minimum virtuality scale Qo — 0.7 GeV 
and the dynamical choice Qo = y/E/L where L is the 
length of the parton path in the medium. As discussed 
in more detail in [l3j . there is no compelling reason to 
identify the Cooper- Frye freeze-out hypersurface in a hy- 
drodynamical model with a sharp surface beyond which 
the medium interaction with partons is absent, and pre- 
scription to let the medium fade out gradually would be 
desirable. However, since the length dependence enters 
under the square root, the dependence on details is suffi- 
ciently weak such that one may hope to get a reasonable 
estimate of the effect. 

As seen from the figure, both computations are able to 
describe Raa well, however the Pt dependence is some- 
what different. This is not surprising, since the expres- 
sion for Qo also contains a dependence on the initial par- 
ton energy E which is proportional to the final hadron 
Pt- Generically, for the dynamical scale choice, a rise of 



Raa with increasing Pt is expected. The curve for fixed 
Qo nas been obtained for a randomized formation time. 
The result computed for average formation time is not 
shown, as it is indistinguishable given the averaging over 
many different in-medium pathlengths as required in the 
computation of Raa- 

Note that the Pt dependence of Raa is different from 
results shown in [l3|- The reason is that in these 
works the matching procedure between hadronic and par- 
tonic scale described in [2!| has not been used, which 
leads to an unrealistic decrease of Raa with Pt- 
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FIG. 4: (Color online) The nuclear suppression factor Raa 
for 20-30% central 200 AGeV Au-Au collisions computed for 
two different options (see text) in YaJEM, compared with 
PHENIX data lldl. Shown are in-plane emission (solid) and 
out of plane emission (dashed). 



The result for 20-30% centrality is shown in Fig. [4l 
separated into in plane and out of plane emission, again 
for a fixed value of Qo and a dynamically computed value 
and compared with PHENIX data fioj . 

It is immediately apparent from the figure that a com- 
putation with fixed Qo results in a spread between in- 
plane and out of plane emission much smaller than seen 
in the data. This is consistent with the observation made 
earlier that the effective pathlcngth dependence of energy 
loss in YaJEM is approximately linear, and the spread 
found agrees with results from models in which the de- 
pendence was assumed to be linear [TJ [l2| ■ 

On the other hand, both spread and normalization of 
the data are well reproduced by a dynamical choice for 
Qo which introduces for light quarks and gluons a strong 
pathlength dependence of energy loss. 

The scenario with dynamically computed Qo has thus 
clear observable consequences when observing hadrons at 
higher energies, such as available at the LHC. Since the 
maximal in-medium pathlcngth is always limited by the 
nuclear diameter times a factor of order one, Eq. ([5]) im- 
plies that Qo will grow with increasing energy, thus the 
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suppression induced by the medium will decrease, leading 
to a rise of Raa with increasing Pp. Somewhat unfor- 
tunately, at LHC kinematics, such a rise is expected in 
various models for a number of reasons, among them the 
shape of the parton spectrum or the fact that the MMFF 
is, even for fixed Qo, energy dependent and reverts for 
large energies to the vacuum result [29| . However, the rise 
induced by the dynamically computed Qo occurs much 
faster than any of these other effects. 

In contrast, the ratio of computed Qo at fixed energy 
for different pathlengths is a constant. When considering 
in-plane and out of plane emission for non-central colli- 
sions, this implies that while the mean (angular averaged) 
Raa changes rapidly with energy, the spread between in- 
plane and out of plane will be more robust and change 
much slower with energy (since asymptotically at high 
energies the MMFF reverts to the vacuum FF, eventu- 
ally the spread must vanish as well). Thus, the scenario 
proposed here can be tested by an observation of in-plane 
vs. out of plane emission of hard hadrons in non-central 
collisions for a sufficiently large range in Pp. Detailed 
predictions for LHC kinematics will be postponed to a 
subsequent publication, as they involve a significant nu- 
merical effort. 



IV. DISCUSSION 

From the results of the previous section, it would ap- 
pear that in the kinematical domain currently experimen- 
tally accessible, imposing a global virtuality resolution 
scale Eq. ([5]) down to which the medium can influence 
the shower development has a much larger influence on 
the observable pathlength dependence of high Pt hadron 
production than a local coherence condition imposed for 
each branching which accounts for LPM interference. 

Note that Qo acts in the presence of a medium quite 
different than in the vacuum case. In vacuum shower, 
Qo regulates the transition from a perturbative partonic 
description in terms of a shower evolution to a non- 



perturbative hadronization description in terms of the 
Lund model, and there is some transition region in which 
both prescriptions lead to similar results such that the fi- 
nal result is not strongly influenced by the detailed scale 
choice. However within YaJEM, by construction (and 
due to a lack of insight into what physics hadronization 
inside the medium would entail), the medium only in- 
fluences the partonic dynamics and there is no compen- 
sation of e.g. a high scale choice leading to weak show- 
ering by an increased action of the Lund model as in 
the vacuum case. Thus, here Qo more directly regulates 
what part of the shower evolution happens in medium 
and what needs to be taken care off by vacuum dynam- 
ics later. 

If this result is confirmed by more detailed investiga- 
tions, it would imply the the idea to use Raa{<P) to dis- 
tinguish between pQCD and strongly coupled models on 
the basis of quadratic vs. cubic pathlength dependence 
is largely a red herring, as the leading effect seen in the 
data is something entirely different. However, since the 
scenario also predicts an energy dependence of the vir- 
tuality scale, translating e.g. into a rise of Raa with 
Pt, testing it should be comparatively easy with a large 
enough kinematic lever-arm as e.g. the LHC experiments 
can provide. 

In general, it seems that in the context of MC shower 
models which can incorporate finite energy effects and to 
some degree also quantum interference, the question of 
how pathlength dependence arises and how to define it 
from obscrvablcs becomes significantly more complicated 
than in the asymptotic limits of analytical models. 
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